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Oxidization of a dilute Si(Gc) alloy is modeled using an original protocol based on molecular dynamics 
simulation and rules for the oxygen insertions. These rules, deduced from ab-initio calculations, favor the 
formation of Si02 against Ge02 oxide which leads to segregation of Ge atoms into the alloy during the 
oxidization front advance. Ge condensation is then observed close to the Si02/Ge interface due to the 
strain induced by oxydization in this region. From the analysis of the simulation process, we propose a one- 
dimensional description of Ge condensation wich perfectly reproduces the evolution of the Ge concentration 
during oxidization of the SiGe alloy. 
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In the challenging world of nanotechnologies, the so- 
called, Ge-condensation technique, is a promising way to 
realize thin Ge layers of good quality for high mobility 
MOSFET channel^ . The method consists in oxidizing at 
high temperature an Sii-jjGe^; alloy epitaxially grown on 
SOI substrate and preliminary capped by a thin silicon 
layer. After the oxidization of the silicon cap that forms 
Si02 oxide, the pursuit of oxidization induces the seg- 
regation of Ge atoms in the SiGe alloy which becomes 
enriched in Ge. 

The development of predicting numerical tools to op- 
timize manufacturing steps^ is essential and necessary 
to get a good understanding of the relevant mechanisms 
involved during the Ge-condensation process. In the con- 
text of the continuous decrease of the microsystem sizes 
and the use of smaller quantities of matter, descriptions 
and precise characterizations at the atomic scale is a 
necessary road map to progress in the understanding of 
the driving mechanisms involved in the Ge condensation 
technique. Such an atomistic description is a challenging 
problem which needs to deal with a complex system and 
many mechanisms occurring at the same time. In par- 
ticular, it requires to model the oxidization of the SiGe 
alloy and the segregation of Ge atoms from the oxydised 
region to the substrate. A full description taking into 
account all mechanisms and reproducing the time evo- 
lution of the system in one model is obviously far from 
the actual state-of-art of computational science at atomic 
scales. 

In this letter, we propose an atomistic model which 
couples molecular dynamics simulations and rules de- 
duced from ab-initio calculations to model oxidization 
of SiGe alloy involving the Ge segregation. In order to 
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be able to describe a system sufficiently large to be rep- 
resentative of the Ge condensation, we use the Stillinger- 
Weber (SW) empirical potential which has been success- 
fully employed and optimized to describe SiGe alloys and 
Si-0 mixed systems by Molecular Dynamics or Monte 
Carlo simulation'^"^. The parameters of this potential, 
which is composed of 2-body and 3-body interactions, 
were fitted to ab-initio calculations and thermodynami- 
cal properties of Si, Ge, and SiGe compounds. For the 
mixed Si-Ge-0 system, we derived a set of parameters 
from those developed by Watanabe et alJ and Clancy ct 
al.^ to study respectively Si-0 and a SiGe. In practice, for 
the Ge-0 interactions, we scaled the SW parametrization 
of Watanabe et al. by using the size and energy param- 
eters of the Ge parametrization developed by Clancy et 
al.^. Arithmetical and geometrical averages have been 
used for the size and energy parameters respectively the 
in 3-body interactions involving triplets composed of Si, 
O and Ge atoms such as OSiGe, SiSiGe, SiOGe, ... 

As mentioned above, the Ge condensation process be- 
gins by the oxidization of the silicon cap covering the 
SiGe alloy. This step is essential since it initiates the 
formation of a Si02 oxide^°. After the oxidation of the 
whole silicon cap, one would expect the formation of ei- 
ther Si02, Ge02 or a mixed Si02-Ge02 amorphous ox- 
ide. However, the experiments show a preferential oxi- 
dization of the Si atoms which induces the segregation 
of Ge atoms and an increase of the fraction of Gc in 
the SiGe alloy^'^^. This preferential oxydization, can be 
understood from the differential in formation energy Ef 
between Si02 and Gc02 in the a-quartz phases : 

Si + O2 ^ Si02 
Ge -t- O2 — > Ge02 

From the total energies obtained using the Wien2k 
code^^ the formation energy can be calculated as : 

Ef°^ = E^°^ -E^ - E°^ (1) 
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where E-^'^^ is the total energy of the oxide, is the 
energy of X = Si or Ge in the diamond phase structure 
and i?*^^ the energy of an isolated O2 molecule. We have 
^ is much lower than E'j'^'^^ , respec- 

8.2 eV and -4.7 eV^'^, in full qualitative agreement 
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with experimental observation since A = E' 
= 3.5 eV implies that the formation of Si02 
favored against that of Ge02- 

The vality of our parametrization of the semi-empirical 
SW potential was also checked on its ability to describe 
this preferential oxydation. The difference between the 
formation energy of Si02 and Ge02 was found to be 
A = 1.8 eV, which is half the value obtained by the 
ab-initio calculations, but widely sufHcient to ensure the 
empirical model to clearly favor the formation of Si02 
against GCO2 phase. 





FIG. 1. (color online) Scheme of the oxygen insertion taking 
into account the preferably oxidization of the Si atom i) the 
highest Si in the substrate surrounded by at least one and 
at most three oxygen atoms is located (yellow circle labeled 
1), ii) when the last chance to introduce oxygen is on a Si- 
Ge bond, the considered Ge atom (bleu circle numebered 2) 
is exchanged with a non-oxidized Si atom of the substrate 
(yellow circle labeled 3). 

The simulation box that we have used was oriented 
in such a way that the x-, y-, and z-axis of the simula- 
tion box are along the <100> directions of the diamond 
lattice. The system size is 27.15 x 27.15 x 108.6 A^. Pe- 
riodic boundary conditions were applied in the [100] and 
[010] directions {i.e. x and y) and the system was free to 
relax in the [001] direction (free surfaces). The starting 
configuration was a bare Sio.95Geo.05 substrate (3700 Si 
and 194 Ge atoms), which was capped by a pure Si layer 
in the upper part (300 Si additional atoms). 

To simulate SiGc oxydization, we use a protocol based 
on molecular dynamics (MD) simulations and on the one- 
by-one incorporation of oxygen atoms into Si-Si bonds. 
This protocole is an extension to the one that we have 
previously used in the case of Si oxydization^. We 
take now into account the preferential formation of Si02 
against Ge02 oxide. In order to simulate the initial nu- 
clei of the oxidization process, oxygen atoms are placed 
on 10 % of the Si dangling bonds (randomly chosen) on 
top of the cap Si layer, which is not reconstructed before 
the oxygen insertion. This initial configuration is equili- 
brated at 1200 K during 20 ps by rescaling the atomic ve- 
locities every 0.5 ps. Then, we include the oxygen atoms 
one by one following the following sequence. We search 



for the highest Si atom in the substrate which is bound 
to at least one and at most three oxygen atoms and the 
oxygen atom is placed in the middle of the longest Si-Si 
bond available around this Si atom^. When only a Si-Ge 
bond is found, we place the oxygen atom in the middle 
of the bond and we exchange the Gc atom with a non- 
oxidized Si atom of the SiGe substrate. The latter being 
determined in order to be the most energetically favor- 
able exchange. This procedure is schematized in the Fig. 
1. The whole system is relaxed by a Molecular Dynamic 
run at 1200 K during 1 ps. The full scheme (insertion, 
relaxation) is repeated for 6000 oxygen atoms insertions 
during which some 600-700 Si-Ge exchanges were real- 
ized. This procedure leads to the formation of a pure 
amorphous Si02 oxide which is composed of a random 
network of Si04 entities connected by vertex, which is 
very similar to that of the Si02/Si system^. 

The Ge concentration profile along the z-axis below 
the Si02/SiGe interface for a typical run is shown in Fig. 
2 after different oxidation steps (5, 2000, 4000, and 6000 
oxygen insertions). The concentration is defined for dif- 
fuse layers of thickness Az = 1 A. Starting from an uni- 
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FIG. 2. (color online) Evolution of the Ge concentration 
?^GG/("■GG + nsi) profile for the Si02/SiGe system after in- 
sertion of 5, 2000, 4000, and 6000 oxygen atoms. O, Si and 
Ge atoms are represented in red, green and black respectively. 
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form concentration of Ge along the z-axis (Fig. 2(a)), one 
can see that the Ge atoms concentrate in the vicinity of 
the interface between Si02 and SiGe (Fig. 2(b)-(d)), in 
good agreement with RBS measurements^^. 

The overal picture is that the oxidation evolves by 
pushing a Ge strip which thickens all along the process, 
whilst preserving a concentration lower than 50 % per 
plane. This means that the Go atoms expelled from the 
oxide mainly exchange with Si substrate atoms located 
close to the Si02/SiGe interface, even though a few ex- 
changes also involve Si atoms far from this interface, lead- 
ing to a slight increase of the overall Ge concentration. 

This qualitative behaviour can be confirmed by ana- 
lyzing the relative position of the Si and Go atoms which 
are exchanged during the oxidation process. In Fig. 3 
we show an histogram of the projections of the distance 
parallel d\\ and perpendicular d± to the interface. The 
distribution of d|| is uniform and does not exceed half 
the size of the simulation box, which indicates that ex- 
changes occur between Si and Gc atoms ranging for all 
the simulation box width. On the contrary, the distribu- 
tion of d± presents a narrow peak which is almost com- 
pletely damped beyond 5 A which corresponds to three 
(001) interatomic layers. This shows that the Ge atoms 
mainly exchange with Si atoms located in the vicinity of 
the Si02/SiGe interface, which leads to their accumula- 
tion (condensation) in this limited region. The tail of the 
distribution of d± towards larger distances confirms that 
a few exchanges can also occur between more distant Si 
and Ge atoms, leading to the slight increase in the overal 
Ge concentration profile. 
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FIG. 3. (color online) Distribution of the parallel dy and 
perpendicular d± projections of the distances separating ex- 
changed Si and Ge atoms during the oxidation process. 

It remains to identify the physical origin of the prefer- 
ence for Si-Ge exchanges to occur close to the interface. 
A carefuU analysis of the local concentrations does not 
allow to establish any correlation between the Ge con- 
centration in the SiGc substrate and the new local envi- 
ronment of the exchanged Ge atom. However, we have 



found that the this preference is induced by the strain 
originated by the oxidation in the interfacial region. In 
Fig. 4 we show the pressure maps averaged over the first 
three SiGe layers below the oxidation front after inser- 
tion of 2000, 4000 and 6000 oxygen atoms. These two 
dimensional maps, parallel to the interface, are very sim- 
ilar at each step of the process and reveal an inhomoge- 
neous strain distribution which is rapidly damped beyond 
these three layers. More precisely, one sees that tensile 
and compressive regions coexist in almost equal propor- 
tions. The Ge condensation close to the interface is then 
originated by the segregation of the " big" Ge atoms (the 
size-mismatch between Si and Ge is of about 4 %) in 
the tensile regions. The number of tensile sites decreases 
when the Ge concentration increases and should disap- 
pear when the concentration reaches 50%. This is indeed 
what is observed in Fig. 2(d), confirming the strain origin 
of Ge condensation close to Si02/SiGe interface. 
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FIG. 4. (color online) Lateral stress maps of the first three 
layers below the Si02/SiGe interface after insertion of 2000, 
4000 and 6000 oxygen atoms. Tensile and compressive regions 
are represnted in blue and red respectivelly. 

From the previous analysis, it is clear that one of the 
essential ingredients which are responsible of the Ge con- 
densation, is the peculiar distribution of distances be- 
tween exchanged Si-Ge atoms displayed in Fig. 3). A 
simple one dimensional model of Ge condensation can 
be used to understand the role played by the over- 
all shape of this distribution. Let us consider a non- 
oxidized SiGe substrate composed of N layers labelled p 
(jj — I, N) , each layer contains a finite number of atoms 
n{p) = nQc{p) + nsiip) where UQcip) and risi(p) are re- 
spectively the number of Gc and Si atoms in the p layer. 
Assuming that oxidation is progressing layer-by-layers, 
by exchanging Ge atoms from the first p layers, starting 
fromp = 1, with Si atoms of deeper layers, distant by d±^ 
from the former, the Ge condensation follows the simple 
set of equations : 



nceip) = nceip) - 1 (2) 

nsiip) = nsiip) + 1 (3) 

riGeip + d\_) ^ nGe{p + di_) + I (4) 

nst{p + d\_) = nsi{p + di_) - I (5) 



Starting from the same initial Ge concentration 5 %, 
number of atoms per layer n{p) = 50 and number of SiGe 
layers N — 72 than in the MD simulation, the Equations 
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(2)- (5) lead to the evolution of the Ge concentration pro- 
files displayed in Fig. 5. The inset in this figure recalls 
the probability distribution of the interlayer distance d± 
between the exchanged Si and Ge atoms which is used 
in the simulation, i.e., for each exchange we randomly 
choose a d± following this distribution. 
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FIG. 5. Ge concentration after each four layers depleted in 
Ge. The inset shows the distribution of SiGe exchanged dur- 
ing the simulation. 



As can be seen, this simple model reproduces the es- 
sential features of the MD model, namely the formation 
of a Ge strip below the Si02/SiGe interface, which thick- 
ens as oxidation goes along. However, the model misses 
the saturation of Ge segregation at about 50 % per plane. 

To summarize, we have shown here by using atomistic 
models combining ab-initio and empirical potentials), 
that the main driving force for Ge condensation under 
oxidation are the preferential Si02 formation against that 
of Gc02 and the existence of an inhomogeneous strain ex- 



tending on a few layers below this interface. The latter 
effect is at the origin of the localisation of the Ge agre- 
gates close to the interface. In addition, we have shown 
that a simple ID model, based on the single knowledge 
of the distribution of distance between Ge and Si atoms 
exchanged during the simulation (related to the strain 
distribution), is sufficient to account for the main fea- 
tures of Ge condensation, providing us a fiexiblc method 
to study this phenomenon under various conditions. 
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